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ABSTRACT 

We present photometric and spectroscopic analysis of the bright detached ecUpsing 
binary BG Ind. The masses of the components are found to be 1.428 ± 0.008 and 
1.293 ± 0.008 Mq and the radii to be 2.290 ± 0.017 and 1.680 ± 0.038 Rq for primary 
and secondary stars, respectively. Spectra- and isochrone-fitting coupled with colour 
indices calibration yield [Fe/H] = -0.2 ± 0.1. At an age of 2.65±0.20 Gyr BG Ind 
is well advanced in the main-sequence evolutionary phase - in fact, its primary is at 
TAMS or just beyond it. Together with three similar systems (BK Peg, BW Aqr and 
GX Gem) it offers an interesting opportunity to test the theoretical description of 
overshooting in the critical mass range 1.2 - 1.5 Mq. 
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1 INTRODUCTION 



The 6th maRnitude star B G Ind was classified as an F3 
dwarf by iMalarodal (Il975l ). an d several years later estab - 
lishe d as an eclipsins binary by [Manfroid fc Mathvd (|l984l ) 



and iMathys. Manfroid fc Renson ( 19861 ). who identified it 
as a detached system with partial eclipses and a period of 
1.464047 d based on ubvy data from 1984. The first ra- 
dial velocity measurements of BG In d were performed by 
I Andersen. Jensen fc NordstromI (| 19841 ). As they wrote, "The 
plates immediately showed spectral lines of two components 
of rather similar type, the lines of one component (the pri- 
mary) being noticeably stronger and broader than those 
of the secondary." Although their data were rather scarce 
(just two spectra), the resulting estimates of masses and 
radii of the components (mi ~ 1.4 Mq, Ri ~ 2.0 Rq and 
m.2 ~ 1.2 Mq, 7?2 ~ 1.5 Rq) proved to be surprisingly 
accurate. Additional Stromgren ph o tometry was collected 
in 1986 by I Van Hamme fc ManfroidI l|l988: hereafter VHMI ) 
who also solved the vby light curves of the system, and found 
an improved value of the period {P — 1.464069 d). For the 
fixed mass ratio q = 0.85 taken from Andersen et al. (1984) 



* Based in part on data obtained at the South African Astro- 
nomical Observatory, 
f E-mail: mnr@camk.edu.pl 
X for the Pi of the Sky team 



they obtained mi = 1.41 Mq, Ri = 2.22 Rq and m2 = 1.20 

Mq, R2 = 1.60 Rq. 

Revised data of iMathvs. Manfroid fc RensonI (|l986l ) 
and VHM, toget her with additional p oint s from 1987, 
were catalogued bv lManfroid et aL I (|l99ll ) and lSterken et al.l 
(1993) as a part of the Long Term Photometry of Vari- 
ables program conducted at ESQ. BG Ind was moni- 
tored by Hipparcos satellite Ferryman et al.l (|l997l ) and 
robotic tel escopes ASAS (iPoimanskil l200ll ) and Fi of 
the Sky (|Malek et all |2010| ): it is also included in 
The Geneva-Copenhagen survey of the solar neighborhood 
l|Holmberg. Nordstrom fc Andersenll2009l ). Recently, a spec- 
troscopic solution of BG In d based on data co llected in 
2006 has been published by iBakis et al.l l|2010t ). who for 
P = 1.464069 d and i = 74.14° found by VHM obtained 
mi = 1.47 ±0.01 Mq, m2 = 1.31 ±0.01 Mq, and a semima- 
jor axis a = 7.64 ± 0.04 _Rq. 

With so much data available, BG Ind may seem to be 
well explored and hence of little interest. However, there 
are at least three reasons to investigate it more thoroughly 
than in the papers mentioned above. First, VHM encoun- 
tered problems with phasing, and they recommended "fur- 
ther monitoring of this binary system for minimum times in 
order to obtain improved ephemeris, and to allow a study of 
the behavior of its orbital period". Second, systemic radial 
velocity measurements gave diverging results: = 39.8 ± 4 
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km s"^ (VHM) and 59.4 ± 5 km s"^ l|Bakis et al.ll20ld ). 
To make the confusion even larger, from three velocity 
components of BG Ind with respect to the Sun listed by 
iHolmberg. Nordstrom fc AndersenI (|2009l ) one obtains the 
total velocity of 20.3 km s~^. 

Third, and most important, both components of BG 
Ind belong to the interesting mass range of 1.1-1.5 Mq in 
which convective cores begin to develop, affecting evolution- 
ary tracks and isochrones via overshooting-related uncer- 
tainties. Moreover, their masses are remarkably similar to 
those of BK Peg - a main sequence binary recently stud- 
ied by [a ausen et al.l (|2010l ). Since the metallicity is similar 
in both cases, and smaller stellar radii observed in BK Peg 
[Ri = 1.987 ±0.008 Rq,R2 = 1.473 ±0.017 Rq) indicate a 
slightly less advanced evolutionary stage, these two systems 
are a potential source of valuable information concerning 
the endphases of hydrogen burning in stars somewhat more 
massive than the Sun. Bearing this in mind, we decided to 
re-analyze BG Ind with the aim to determine its precise 
physical parameters and evolutionary status. 

Our paper is based on photometric and spectroscopic 
data described in Sect. 2. The analysis of the data is detailed 
in Sect. 3, and its results are discussed in Sect. 4. 



2 OBSERVATIONAL MATERIAL AND DATA 
REDUCTION 

2.1 Photometry 

Our photometric solutions are based on six sets of photo- 
metric measurements listed in Table [T] spanning a period 
from May 1986 to November 2009 (points with Am > 0.1 
mag, where Am stands for the deviation from the prelimi- 
nary fit (see Sect. 13. 2| . were rejected from the original ASAS 
and Pi datasets). We found that all photometric data from 
Table [T] neatly phase with the ephemeris 



to = HJD 47876.3792 ± 0.0004 
P = 1.46406335 ± 0.00000002 d 



(1) 



(see Fig. IA1|I , and we are rather confident that the period 
of BG Ind remained constant for over 23 years. The only 
indication suggesting a possible period change comes from 
the earliest observations of this system (runs 1 and 2 of 
VHM, not included in our analysis). 

Available online are als o Hipparcos Bt and Vt lightcurves 
llPerrvman et all Il997l ) , and Stromgren u l ightcurve of 
iManfroid et al.l (|l99ll ) and ISterken et all (|l993l ) . Their qual- 
ity was too poor to use them for photometric solutions; how- 
ever Bt and Vt data yielded a useful estimate of the tem- 
perature of the hotter c omponent (see Sect.l3.2[>. We note 
here that according to (|Suchkov. Makarov fc VogesI |2003| ) 
the Stromgren (6 — y) excess of BG Ind amounts to 0.001 
only, so that reddening effects can be neglected. 



2.2 Spectroscopy 

Because of problems with the systemic velocity of BG Ind 
mentioned in Sect. [1] we decided to use our own spectral 
data, collected in September / October 2007 with the fiber- 
fed Giraffe spectrograph on the 1.9-m Radcliffe telescope 
at the South African Astronomical Observatory. The seeing 



Table 1 


List of light 


curves used 


in this paper for photometric 


solutions of BG Ind. 










Label 


Filter 


1=* day 


last day 


Number of 


Ref. 






[HJD-2400000] 


data points 




MSb 


Stromgren b 


44582 


47070 


165 


1,2 


MSv 


Stromgren v 


46582 


47070 


172 


1,2 


MSy 


Stromgren y 


46582 


47070 


172 


1,2 


HipH 


Hipparcos H 


47888 


49054 


142 


3 


ASV 


Johnson V 


51997 


55167 


576 


4 


PiR 


Cousins R 


53902 


54933 


828 


5 



3: 



1: IManfroid et al.l l ll99ll'l: 2: 'Sterken et al.' ('l9 9j; 
IPerrvman et al.l lll997h : 4: IPoimansk i (.2001); 5: Malek et al.l 
||2010|'I . Pi of the Sky observations were done wihtout filters, but 
the combined sensitivity of optical system and detector closely 
matched that of Cousins R-filter. 



oscillated between 1.5 arcsec and 2.5 arcsec. A 2.7 arcsec 
entrance window provided a resolution of almost R = 40000 
at A = 4470 A. During the observations pairs of scientific 
spectra were taken, separated by an exposure of a thorium- 
argon hollow-cathode lamp. The exposure times per spec- 
trum ranged from 400 s to 900 s. The observations were 
reduced within the IRAF Q ECHELLE package. After bias 
and fiat-field correction each pair of the frames was com- 
bined into a single frame, allowing for the rejection of cos- 
mic ray hits. Altogether, 23 reduced spectra were obtained. 
For further analysis a wavelength range extending from 
4300A to 5800A was used, in which most of the spectra 
had 20 < S/N < 40 with a few cases of lower quality. 

Our method of radial velocity measurements is based on 
the broadening function formalism introduc ed by Rucinskil 
( 2002) and it was described in detail by iKaluznv et al.l 
(|2006l i.FlThe spectrum of the bright star HD 200163 (spec- 



tral type F3V), obtained in the same observing period and 
with the same instrumental setup, served as a template. 
The resulting barycentric radial velocities are listed in Table 
[21 The mean error, estimated from velocity-curve fitting, is 
±0.63 km s-i 

The binary sp ectrum of BG Ind was disentangled using 
the code of Konac ki et al.l ()2010l ) with the aim to estimate 
the temperature of the components. Unfortunately, because 
of strong rotational broadening and blending of the lines 
the only method we could apply was a direct comparison to 
synthetic spectra (see Sect. 13.2]) . 



3 ANALYSIS 

The compo nents of BG Ind h ave masses in the range 1.1-1.5 
Mq fVHM. ISakis et al.ll2010l '). so that their envelopes should 
be convective. Consequently, we adopted gravity darken- 
ing exponents gl = g2 = 0.32 and bolometric albedos 
Al = A2 = 0.5. Just in case, we also obtained a solution 
with radiative envelopes (gl — g2 — Al — A2 — 1.0) and 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the AURA, Inc., under coopera- 
tive agreement with the NSF. 
^ The software used in this paper is 



this paper 
|http://users. camk.edu.pl/pych/BF/ 1 



available at 
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Table 2. Barycentric radial velocities of BG Ind phased according 



to the ephemeris 


given 


3y equation ([Tjl. 




HJD-2454000 


Phase 


fi [km s~-^] 


V2 [km s ■'■] 


363 


398889 


U-OOO 


134 234 


-00. ( 0\J 


363 


420815 




1 97 71 9 


7fi i^l 


366 


290933 


U.OUO 


142 702 


QI A'^P. 
-y 1 .^ou 


368 


315269 




- / o.oy 1 


1 1^9 7A(^ 


368 


335724 




1 7^ 


1 1;7 /iQ/i 


368 


365740 




07 91^9 


iOi .OiU 


369 


244426 




-LOU. 000 


91 )R 


369 


263211 


u.ooo 


101 . 1 Ui 


7Q fl'^7 

- ( y .00 / 


371 


249086 




Qfl (\^(\ 
-OU .\jO\j 




371 


266822 


n 907 
U-ZU ( 


QO 99/1 


1 ^7 009 


371 


284471 


n 91 Q 


1 Q7 


1 f^Q f^flQ 

ioy.ouy 


371 


315563 


0.240 


-86.278 


161.496 


371 


328299 


0.249 


-86.007 


162.845 


371 


344954 


0.260 


-86.568 


162.832 


371 


361949 


0.272 


-85.781 


160.540 


371 


379009 


0.283 


-83.505 


158.200 


371 


395454 


0.295 


-82.164 


157.302 


371 


412749 


0.306 


-79.292 


153.182 


372 


255302 


0.882 


112.450 


-57.583 


372 


275909 


0.896 


104.275 


-48.576 


372 


292410 


0.907 


96.605 


-41.796 


372 


308893 


0.918 


89.945 


-33.640 


376 


386726 


0.221 


144.054 


-92.932 



we found that even within unrealistically small formal error 
limits calculated by PHOEBE which are all "radiative" pa- 
rameters but one did not differ from the "convective" ones. 
The exception was the radius of the primary Ri which in 
the radiative case was by 5.6% smaller (we define primary 
as the more massive component). 

Since BG Ind is a short-period system, and broadening- 
function fitting yielded v sin i values close to those expected 
from the synchronous rotation for the radii obtained by 
VHM, we assumed full synchronization of both components. 
The efTects of reflection were incl uded, and a l ogarit hmic 
limb-darkening based on tables bv lVan Hammi (|l993h was 
used as impleme nted in PHOEBE 031a. 

According to lHolmberg. Nordstrom" fc AndersenI (|2009l ') 
the system is slightly underabundant in metals {[Fe/H] — 
—0.3). However, their estimate is based on Stromgren pho- 
tometry of the total light, and as such it cannot be 100% re- 
liable. As for the temperatures of the components, VHM as- 
sumed Ti — 7000 K for the primary, and obtained T2 — 6450 
K for the secondary from the fit. While these values are cer- 
tainly reasonable, they may be too high for the radii ob- 
tained by VHM which indicate that at least the primary is 
about to leave the main sequence or even has left it. More- 
over, upon comparing our phased velocity curve with phased 
light curves from Table [T] we found that the hotter compo- 
nent of BG Ind is the secondary (see Fig.[T]).[f| The s ame con- 
clusio n follows from the velocities measured by Bakis et al.l 
(|2010l ) which also nicely phase with the ephemeris IT} , show- 
ing, however, a much larger scatter (the rms residual from 



^ The misidentification of VHM is entirely understandable when 
one remembers that they had just two velocity measurements and 
only an approximate value of the period. 




-a - 



-0.5 0.5 

Figure 1. Velocity curve of BG Ind based on data contained in 
Table[2] Systemic velocity and rms residuum are equal to 31.2 ± 
0.1 km s~^ and 0.63 km s~^ respectively. Phase corresponds to 
the center of the shallower photometric minimum. 

our fit to their data is 4.6 km s~^, only slightly larger than 
4.4 km s~^ found in the original paper). 
Thus, while searching for an observational estimate of the 
temperature of BG Ind components we decided to focus on 
the secondary. As it rotates significantly slower than the 
primary and must be less evolutionary advanced, we ex- 
pected it to be a more or less normal main-sequence star 
to which available color-temperature calibrations could be 
reliably applied. 

3.1 Preliminary solutions 

We started the analysis from searching for solutions with T2 
equal to 7000, 6500 or 6000 K, an d \Fe/H] equal to 0.0 or 
-0.3. We used PHOEBE in terface llPrsa fc Zwitt"eill2005l ') to 
the Wilson-Devinney code jWilson fc Devinnevlll97ll ). solv- 
ing for all lightcurves simultaneously. The essential aim of 
preliminary calculations was to check how strongly in this 
parameter range vary the gravitational acceleration of the 
secondary g2 and the contributions of the secondary to the 
total light in various wavebands at ip = 0.25. As detailed in 
the next Section, the first of those parameters was needed 
for an estimate of T2 based on our spectra, while the re- 
maining ones - for an independent estimate of T2 based on 
calibrations of color indexes. 

There is no indication for a nonzero eccentricity in ei- 
ther light or velocity curves, but just in case we included 
iterations of e. The solutions yielded e — 0.0 ± 0.0005, and 
an almost constant log 32 = 4.14 — 4.15. The primary's accel- 
eration gi was also almost constant, but smaller by a factor 
of 2, i.e. appropriate for the beginning of the subgiant branch 
rather than for the main sequence. 

3.2 The temperature of the secondary 

With 32 fixed, we generated an array of synthetic spectra 
for 6000K s; T ^ 7000K and - 1 < [Fe/H] < based on the 
library of (|Coelho et al.|[2005l l. The respective spacings AT 
and A[Fe/H] were equal to 100 K and 0.1. All spectra were 
rotationally broadened with vsini = 78 km s^^ for the pri- 
mary and 53 km s"'^ for the secondary, consistently with the 
mean values of i and component radii obtained from prelim- 
inary fits (73°, 2.35 Rq and 1.61 Rq, respectively). Next, 
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Table 3. PHOEBE solutions with T2 = 6650 K and [Fe/H] = -0.2. 



light 


rms 


i 


a 


mi 


7712 


Ri 


R2 


Ti 


Li 


L2 


curve 


[mmag] 


[deg] 


[Rq] 


[Me] 


[Mq] 


[Re] 


[Re] 


[K] 


[Lq] 


[Lq] 


MSb 


7.0 


72.95 


7.567 


1.428 


1.293 


2.313 


1.651 


6312 


7.37 


4.62 


MSv 


7.5 


72.95 


7.567 


1.428 


1.293 


2.303 


1.663 


6320 


7.34 


4.69 


MSy 


7.9 


73.00 


7.565 


1.427 


1.292 


2.165 


1.796 


6443 


7.00 


5.47 


HipH 


6.4 


73.12 


7.560 


1.424 


1.289 


2.226 


1.723 


6382 


7.13 


5.04 


ASV 


16.1 


72.72 


7.577 


1.434 


1.298 


2.395 


1.637 


6375 


8.22 


4.55 


PiR 


11.4 


73.03 


7.567 


1.428 


1.293 


2.340 


1.610 


6262 


7.30 


4.39 



for each synthetic spectrum the sum of squared deviations 
from its observed counterpart was calculated. In the case 
of the secondary the smallest sum was obtained for T2 = 
6500 K and [Fe/H] = -0.2 (see Fig.\M\ioT the comparison 
of best-fitting spectrum to the observed one). In the case of 
the primary no unique minimum was found, for what broad 
and/or blended lines are probably to be blamed (e.g. Mgl 
lines at 5167A and 5173A were observed as a single spectral 
feature). We estimate the uncertainties of temperature and 
metallicity at ±100 K and ±0.1, respectively. 

In all preliminary solutions with [Fe/H] — the sec- 
ondary to primary luminosity ratio L2/L1 (found from the 
contribution of each component to the total light of the 
system at ip — 0.25) was almost constant in each of b, v 
and y bands, amounting respectively to 0.570-0.573, 0.598- 
0.606 and 0.551-553. Given the total apparent magnitudes 
of BG Ind at y = 0.25, we used these ratios to calculate the 
corresponding apparent magnitudes of the secondary, and 
we found its temperature from (b — y) — T^ff calibration 
of ICasagrande et al.l HoT^). The result for [Fe/H] = 0.0 
was a remarkably constant T2 = 6523 — 6528 K. Solutions 
with [Fe/H] = —0.3 were only slightly more diverging, with 
T2 = 6477-6487 K. After averaging, we got T2 = 6525 ±63 
K and T2 — 6483 ± 63 K, respectively, where the error 
includes inaccuracies of color-temperature calibration and 
magnitu de measurement at (/J = 0.25. Based on th e calibra- 
tions of iHolmberg. Nordstrom fc AndersenI (|2007l 'l. we also 
derived the observed [Fe/H] index, which varied between 
-0.4 and -0.7. While being only marginally consistent with 
the lower of the [Fe/H] values assumed in preliminary fit- 
ting, it confirmed that BG Ind is indeed metal-deficient. 

The third estimate of T2 was based on the Hipparcos 
photometry. Proceeding as before, we found the light ra- 
tio L2/L1 at ip — 0.25 in B t and Vt bands, and, u sing 
(B — V)t — Teff calibration of ICasagrande et al.l (|201Cll ). we 
derived T2 = 6650 ±83 K for [Fe/H] = and T2 = 6496 ±83 
K for [Fe/H] = -0.3. 

For the fourth temperature estimate we used Hipparcos 
parallax of BG Ind (tt = 15.04 ± 1.1 mas) and ASAS pho- 
tometry in l/-band. Proceeding as before and u sing bolo- 
metric corrections of iLeieune. Cuisinier fc Busd l|l998r ). we 
found the observed 1^- band magnitude of the secondary 
(V2°=7.22), and compared it to apparent magnitudes V2 
calculated from preliminary solutions. Solutions with T2 — 
6000 which yielded 7.59 < V2 sj 7.88 were clearly ruled 
out. Those with T2 = 6500 K and T2 = 7000 K produced, 
respectively, 7.27 sC V2 «: 7.56 and 6.95 < V2 < 7.24, i.e. 
magnitudes almost consistent with V2, but remaining on the 
low side in the first case and on the high side in the second 



case. A nearly ideal agreement between observed and calcu- 
lated magnitude was reached in an additional preliminary 
fit with T2 = 6750 K, for which the difference ¥2° - V2 was 
equal to just 0.02 mag. 

Thus, the revised Hipparcos parallax pulled the tem- 
perature estimate firmly upward. Trusting it, we adopted 
T2 = 6650 K with a fiducial uncertainty range of ±100 K, 
marginally compatible with all available data. As it was de- 
tailed above, the trend observed in preliminary solutions in- 
dicated that the metallicities obtained from Stromgren pho- 
tometry would cause T2 to fall too low to be compatible with 
V2 ■ Bearing this in mind, we set [Fe/H] to -0.2, i.e. to the 
value resulting from the spectroscopic estimate. 



3.3 Final model fitting 

For the final analysis the JKTEBOP code (lEtzell [l981 



Popper fc Etzellll98ll:[Southworth. Maxted fc Smallevll2004l : 



Soutworth et al.l 220J) and a spectroscopic data solver writ- 
ten and kindly provided by G. Torres were used except 
PHOEBE. While JKTEBOP gives realistic estimates of the 
errors, it cannot be used when the stars are too distorted. 
Unfortunately, this is the case of BG Ind, where the pri- 
mary's distortion exceeds the allowable limit by 10%. Con- 
sequently, we employed JKTEBOP solely to estimate the 
errors of the photometric solution (note that differences be- 
tween the errors of the slightly inaccurate JKTEBOP solu- 
tion and the errors of the appropriate solution must be small 
quantities of the second order which can be neglected). Simi- 
larly, realistic errors of the spectroscopic solution were found 
using the Torres code. (Because the Torres code requires 
center of mass velocities on input, it was necessary to cor- 
rect the observed light center velocities for effects caused by 
the distortion of the components. Additive phase- dependent 
corrections were calculated using the Wilson-Deviney code; 
their values rang ed from -0.86 to 0.27 km s"^). 

JKTEBOP cannot deal with multiple lightcurves, so 
that the final model fitting had to be performed separately 
for each of the photometric datasets listed in Table [1] We 
fixed e = as indicated by preliminary solutions, and started 
the final calculations from PHOEBE fits which produced 
six sets of parameters listed in Table [3] whose second col- 
umn shows the rms deviation of observed points from the 
fit. Next, we calculated errors of inclination i and relative 
radii of the components ri,2 ~ R\,2/A using JKTEBOP, 
and errors of mi.2sin^i and j4sini using the Torres code. 
The erros were then transformed into errors of parameters 
returned by PHOEBE, and assigned to respective PHOEBE 
solutions. Finally, weighted averages of PHOEBE parame- 
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Figure 2. Location of BG Ind components on M — R, M — logL and logT^jj — logL planes. The lines are Dartmouth [Fe/H] = —0.2 
isochroncs for 2.60 and 2.67 Gyr (i.e. for extreme ages still fitting the primary on the M — R plane). 



ters and the errors of those averages were found from the 
standard formulae 



Table 4. Physical parameters of BG Ind. 



< X > = 



and 



(2) 



(3) 



E:.r 

The final parameters with their errors are listed in Table ID 



4 DISCUSSION AND CONCLUSIONS 

In Sect. [1] we outlined three problems which prompted us to 
analyze BG Ind: variability of the period, doubtful systemic 
velocity and lack of accurate parameters of the system. As 
mentioned in Sect. 12. Il we found that between 1986 and 2009 
the period remained constant at 1.46406335 ± 0.00000002 
d. The only indication that it might have changed comes 
from the earliest observations of BG IND from 1981 and 
1984 (runs 1 and 2 of VHM) which were not included in our 
data. As for the systemic velocity, our value of 31.2 ± 0.2 
km s~^ does not agree with any of those obtained by other 
authors. The value found by VHM (39.8 ± 4 km s"^) was 
based on two measurements only, a nd its error was lik ely 
underestimated. 59.4 ± 5 km s~^ of iBaki^ et al.l (|2010l ) is 
rather large for an F-type star in Sun's vicinity. The origin 
of such a large discrepancy with all remaining estimates is 
difficult to explain - we may only note that they did not 
observe any radial velocity standards, and it is conceivable 
that they reversed the sign of the heliocentric correction 
while reducing the data. Finally, the low Vtotai = 20.3 km 
s~^ of lHolmberg. Nordstrom fc AndersenI l|2009l ) can be ex- 
plained by the fact that they did not detect the binarity of 
BG Ind (this system is not included in their Table 2), and 
their velocity must have been contaminated by the orbital 
motion. 

The parameters found in Sect. 13.31 are accurate enough 
for isochrone fitting. We used solar-scaled Dartmouth 



Parameter 




Value 




Unit 


P 


1.46406335 


± 


0.00000002 


[d] 


i 


72.96 


± 


0.23 


[deg] 


e 


0.000 


± 


0.0005 




a 


7.567 


± 


0.013 


Re 


mi 


1.428 


± 


0.008 


Mq 


m2 


1.293 


± 


0.008 


Mq 


Ri 


2.290 


± 


0.017 


Rq 


R2 


1.680 


± 


0.038 


Re 


logTi[K] 


3.803 


± 


0.018 




logT2[K] 


3.823 


± 


0.015 




logLi/L0 


0.87 


± 


0.07 




logL2/LQ 


0.68 


± 


0.08 




[Fe/H] 


-0.2 


± 


0.1 




age 


2.65 


± 


0.20 


Gyr 



isochrones (Dotter et al. 2008) which include core overshoot- 
ing defined as a product of the pressure scale height and 
a factor Oover which depends on stellar mass and compo- 
sition (at nearly-solar metallicities it grows from 0.05 for 
1.2 s; M s; 1.3M0 through 0.1 for 1.3 < M ^ 1.4Mq to 0.2 
for M > I.4M0). The convection itself is treated according 
to the standard mixing length theory with solar-calibrated 
mixing length parameter a^ni = 1.938. 

What makes BG Ind particularly interesting is that 
the masses of the components fall at the beginning and 
at the end of the range where Oover ramps up. First, we 
fitted isochrones calculated for —0.4 ^ [Fe/H] ^ with 
a step of 0.1 to the most accurately determined parame- 
ters, i.e. masses and radii of the components. The respec- 
tive ages we found were 2.27-2.43, 2.45-2.55, 2.60-2.67, 2.77- 
2.85, and 2.96-3.10 Gyr. Next, we checked how well these 
isochrones perform on M — log L and log Tef / — log L planes. 
The best agreement was obtained for [Fe/H] = —0.2. The 
fit with [Fe/H] — —0.1 was almost equally good; that 
with [Fe/H] = —0.3 was still acceptable; and those with 
[Fe/H] = and [Fe/H] = -0.4 had to be rejected. Thus, 
the isochrone fitting confirmed our spectroscopic estimate of 
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1.2 M/M, 



J \ \ I \ \ I l\ \ I \ \ l_ 

3.85 3.8 3.75 logT^jj 



Figure 3. Systems with similar evolutionary advancement as BG Ind. The lines are Dartmouth [Fe/H] = 0.0 isochrones for 2.60 and 
2.80 Gyr (BK Peg); 2.10 and 2.45 Gyr (BW Aqr); 2.45 and 2.70 Gyr (GX Gem). 



metallicity which we finally fixed at [Fe/H] = —0.2 ± 0.1. 
The corresponding age is 2.65 ± 0.20 Gyr. 

The location of BG Ind on M - R, M ~ logL and 
log Tef f — log L planes is shown in Fig. [2] together with 
t = 2.60 Gyr and t — 2.67 isochrones obtained for [Fe/H] = 
—0.2. One can see that the more massive primary has al- 
most reached the beginning of the subgiant branch, while 
the secondary is still on its way to TAMS. The agreement 
between theoretical and observational data would be ideal 
if it were not for small discrepancies in R2 and Ti (by ~3% 
and ~1.5%, respectively). The first one could originate from 
the fact that the eclipses of BG Ind are partial and the 
secondary is by almost 40% smaller than the primary (and 
therefore less deformed). As a result, and because of rather 



poor quality of available photometric data, the accuracy of 
R2 determination has to be markedly lower than that of Ri . 
The second discrepancy, at a first glance rather insignifi- 
cant, turned our attention because it occurs precisely where 
the effects of overshoot-treatment should be largest (the pri- 
mary of BG Ind is a star with M > IAMq at TAMS). We 
decided to check if the same effect appears in other systems 
with similar masses and in similar evolutionary phase. 

Based on a recent compilation of Clausen et al.l (|2010l ). 
we chose BW Aqr, BK Peg, and GX Gem whose component 
masses range from 1.26 to 1.49 A/0, and whose [Fe/H] in- 
dexes are consistent with 0. The fitting of solar-scaled Dart- 
mouth isochrones for [Fe/H] = yielded respective ages of 
2.15 - 2.45, 2.6 - 2.8 and 2.50 - 2.75 Gyr, placed roughly 
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halfway between Yonsei - Yale a nd VRSS ages quoted in Ta- 
ble 12 of IClausen et all l|2010l '). Figs. [2] and [3] demonstrate 
that the temperature discrepancy, absent in the relatively 
unevolved secondary components of BG Ind and BK Peg, in- 
creases with the evolutionary advancement (whose best indi- 
cator is the distance from the sharp upturn of the isochrones 
on the M — logL plane) until it becomes clearly visible in 
GX Gem whose both components are at TAMS or have al- 
ready left the main sequence. Note that luminosity errors 
are larger than those quoted bv lClausen et al.l (|201(]| ) - this 
is because we recalculated them according to the formula 

to make them consistent with ours. 

The discrepancy is a 1-a effect and as such it may not be 
real, however the trend it exhibits suggests there might be 
some physics behind. It is beyond the scope of our paper to 
identify physical factors or assumptions potentially respon- 
sible for such effect. Whether anybody decides to look for 
them or not, improving the quality of photometric and spec- 
troscopic solutions of all four systems is certainly a worth- 
while task, although in the case of BG Ind it might prove 
rather difficult because of strong rotational broadening. 
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APPENDIX 



Amag 




Figure Al. BG Ind light curves used in tiiis paper (see Table [T] 
for the list), phased with the ephemeris {TJ. Individual curves are 
normalized to magnitude at maximum light. 




Figure A2. A section of the disentangled spectrum of the sec- 
ondary (line) compared with the best fitting synthetic spectrum, 
obtained for T^ff = 6500 K, g = 4.15, [Fe/H] = -0.2, and rota- 
tionally broadened with tisini = 53 km s~^ (dots). The largest 
differences appear in deep lines and in overlap regions of cchelle 
orders. 



